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^ ■ I. INTRODUCTION 

In many physical systems, cooperative interactions between spin-like (two-state) degrees of freedom tend to establish 
some kind of order in the system, while the presence of some noise effect (due to temperature, external transverse 
field etc.) tends to destroy it. Tranverse Ising model can quite succeessfully be employed to study the order-disorder 
\ transitions in many of such systems. 

An example of the above is the study of ferro-electric ordering in Pottasium Dihydrogen Phosphate (KDP) type 
systems (see, e.g., |l|). To understand such ordering, the basic structure can be viewed as a lattice, where in each 
03 ■ lattice point there is a double-well potential created by an oxyzen atom and the hydrogen or proton resides within it 
in any of the two wells. In the corrosponding Ising (or pseudo-spin) picture the state of a double-well with a proton 
at the left-well and that with one at the right-well are represented by, say, | f) and | J.) respectively (see, for a portion 
of the lattice, Fig. 1.1). The protons at neighbouring sites have mutual dipolar repulsions. Hence had proton been a 
classical particle, the zero -temperature configuration of the system would be one with either all the protons residing 
at their respective left-well or all residing at the right-well (corrosponding to the all-up or all-down configuration 
i— 1 1 of the spin system in presence of cooperative interaction alone, at zero-temperature). Considering no fluctuation at 

■ zero temperature, the Hamiltonian for the system in the corrosponding pseudo-spin picture will just be identical to 
the classical Ising Hamiltonian (without any transverse term). However, proton being a quantum particle, there is 

I , always a finite probability for it to tunnel through the finite barrier between two wells even at zero-temperature due 
' to quantum fluctuations. To formulate the term for the tunnelling in the corrosponding spin-picture, we notice that 
\Q ', a x is the right operater. This is because 

o*|T) = U) and ^|4> = |T>, (1) 

m . 

where | |) represents the state where the proton is in the left well, while | J.) represents that with the proton in the 
right well. Hence the tunelling term will exactly be represented by the tranvere field term in the transverse Ising 
Hamiltonian. Here the transverse field coefficient T will represent the tunnelling integral, which depends on the 
width and height of the barrier, mass of the particle, etc. 

a 

^ ; II. TRANSVERSE ISING MODEL (TIM) 



X 



Such a system as discussed above, can be represented by a quantum Ising system, having Hamiltonian 



w = -£ J *K^- r E<- (2) 

Here, is the coupling between the spins at sites i and j, where cr a, s (a = x, y, z) are the Pauli spins satisfying the 
commutation relations 

k")of] = 1i5ije aM a2 (3) 
Here, 5ij is the Kronecker's 5, and e a p-y is the Levi-Civita symbol, and in (1) represents neighbouring pairs. 

The Pauli spin martices being representatives of spin-1/2, a z has got two eigenvalues (±1) corrosponding to 
spins aligned either along z-direction or along the opposite direction respectively. The eigenstate corrosponding to 
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FIG. 1: The double wells at each site (e.g., provided by oxygen in KDP) provide two (low-lying) states of the proton (shown 
by each double well) indicated by the Ising states | f) and | J.) at each site. The tunnelling between the states are induced by 
the transverse field term (Ta x ). The dipole-dipole interaction Jij here for the (asymmetric) choice of one or the other well at 
each site induces the 'exchange' interaction as shown. 



eigenvalue (+1) is symolically denoted by | |), while that corresponding to (—1) is denoted by | |) . 
If we represent 

and 



II) * ( J ) , (4) 
then taking these two eigen-vectors as basis, Pauli spins have following matrix representations 

•■-(?*)■- (f^V-U- .)- < 5 » 

With these, one can see that relations in (3) are easily satisfied and the tunnelling required in (1) can be easily 
accommodated. The order parameter for such a system is generally taken to be the expeectation value of z-component 
of the spin, i.e. (<r z ). Needless to say that in such a system absolute ordering (complete alinement along z-direction ) 
is not possible even at zero-temperature, i.e., (ct z )t=o 1, when r =/= 0. In general, therefore, the order ((a z ) =/= 0) 
to disorder (<j z ) = transition can be brought about by tuning either of, or both of the tunnelling field T and the 
temperature T (see Fig. 1.2). 



III. MEAN FIELD THEORY (MFT) 



(a) For T = 
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Let, 

of = \a\ cos#, and of = |<?| sin#, (6) 

where 9 is the angle between a and z-axis. This renders the two mutually non-commuting part of the Hamiltonian 
(2) commuting, since both are expressed in terms of \a\ operator only. If a is the eigen-value of \a\ (a = 1 for Pauli 
spin) , then the energy per site of the semi-classical system is given by 

E = -aT sin6 -a 2 J(0) cos 2 6, (7) 

J(0) = Ji(0) = Jiji where j indicates the j-th nearest neighbour of the i-th. site. And the average of the 

spin-components are given by 

(a z ) = cos (9 
(a x ) = sm9. 

The energy (7) is minimized for 

sin(9 = ryj(O) or, cos(9 = 0. (8) 

Thus we see that if F = 0, {a x ) = and the order parameter (a z ) = 1, indicating perfect order. 

On the other hand, if T < J(0), then the ground state is partially polarized, since none of (er 2 ) or (a x ) is zero. 
However, if T > J(0), then we must have cos# = for the ground state energy, which means (a z ) — 0, i.e., the state 
is a completely disordered one. Thus, as T increases from to J(0), the system undergoes a transition from ordered 
(ferro)- phase with order parameter (cr z ) = 1 to disordered (para)-phase with order parameter (a z ) — (see Fig. 1.2). 



(b) For T ^ 

The mean field method can also be extended to0, 0] obtain the behaviour of this model at non-zero temperature. 
In this case we define a mean field hi at each site i, which is, in some sense, a resultant of the average cooperative 
enforcement in z-direction and the applied transverse field in x-direction. Precisely, we take, for general random case, 



Yx+\-yjM)\z, (9) 
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and the spin- vector at the i-th site follows h^. The spin- vector at i-th site is given by 

Si = afx + cr-z, 

and Hamiltonian thus reads 

H = -J2^- ( 10 ) 

i 

For non-random case, all the sites have identical ambience, hence hi is replaced by h = Fx + (<r z )J(0). And the 
resulting Hamiltonian takes the form 

i 

The spontaneous magnetization can readily be written down as 

CT = tanh(/3|h|).i- 
|h| 

|h| = v/r2 + (j( )(^)) 2 . (it) 

Now if h makes an angle 9 with z-axis, then cos# = J(0)(a z )/|h| and s'm9 — T|h|, and hence we have 

= \h\co S e = [tanh(/3|h|)] (^p) > 

and 

= [tanh(/3|h|)]^. (12) 
|h| 

Here, j3 = (l/fc^T). Equation (12) is the self-consistency equation which can be solved or graphically or otherwise, to 
obtain the order parameter (<r z ) at any temperature T and transverse field T . Clearly, the order-disorder transition 
is tuned both by T and T (see Fig. 1.2). 

r = (Transition driven by T): 

Here, 

and 

(a*) = 

One can easily see graphically, that the above equations has a nontrivial solution only if fc^T < J(0), i.e., 

(a z ) ^ for k B T < J(0) 

(a z ) = for k B T > J(0). 
This shows that there is a critical temperature T c = J(0) above which, there is no order. 

For k B T — > (Transition driven by T): 
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Here, 

.J(C\)(rr z ) / ^ , 



(Q= - J(0)(f72) (since, tanh, 



a;— »oo 



From this equation we easily see that in the limit T/J(0) — > 1, the only real nontrivial solution is 

- o 

and 

r r 

(a x ) = - -» 1, as — - -» 1. 

V(r) 2 + (J(0)(a-))2 J(0) 

Thus we see that their is a critical transverse field T c — J(0) such that for any r > T c there is no order even at 
zero temperature. In general one sees that at any temperature T < T c , there exist some transverse field T c at which 
the transition from the ordered state {(<r z ) 0) to the disordered state ((o z ) = 0) occurs. The equation for the phase 
boundary in the (T — T) - plane is obtained by putting (a z ) — > in equation (12). The equation gives the relation 
between r c and T c as follows 

'- h G&) = w < i3 » 

One may note that for ordered phase, since (a z ) ^ 0, 

tanh(/3|h|) = —7— = Constant. 
|h| ' J(P) 

Hence, (a x ) = (r/\h\\ tanh(,9|h|) = T/J(0); independent of temperature in the ordered phase. While for the 
disordered phase, since (<j z ) = 0, 

(a x ) = tanh(/3r). 

Using magnetic mapping, mean field theory of this type was indeed applied to (the BCS theory of) superconductivity 
0, as shown in appendix A. 

IV. DYNAMIC MODE-SOFTENING PICTURE 

The elementary excitations in such a system as described above are known as spin waves, and they can be studied 
using Heisenberg equation of motion for a z using the Hamiltonian. The equation of motion is then given by 

of = (i»)- 1 [of,W] (14) 

or, 

of = 2Va\ (with h = 1) 

Hence, 

of = 2Tof = 4r J2 Jv°t°T - 4r 2 of . (15) 

3 

With Fourier transforms and random phase approximation (of o| = of (o|) + (of )o| , with (a z ) = in para phase), 
we get 

^ 2 = 4r(r- J(q)(* x )), (16) 

for the elementary excitations (where J{q) is the Fourier transform of Jij ). The mode corresponding to (q = 0) 
softens, i.e., loq vanishes at the same phase boundary given by equation (13). 
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FIG. 3: The Suzuki- Trotter equivalence of quantum one dimensional chain and a (1+1) dimensional classical system. J' 
indicates the additional interaction in the Trotter direction. 



V. SUZUKI-TROTTER FORMALISM 



Exact analysis for the quantum fluctuation can indeed be tackled by using, renormalization group theory; see 
appendix B for real space quantum RG theory for one dimensional chain (cf However, such formalisms have 
serious limitations in applicability and the Suzuki- Trotter formalism to map the quantum problem to a classical one 
has been of enormous practical importance (e.g. in simulations). 

Suzuki- Trotter formalism is essentially a method to transform a ci-dimensional quantum Hamiltonian into a 
(d+l)-dimensional effective classical Hamiltonian giving the same canonical partition function. Let us illustrate this 
by applying it to transverse Ising system. We start with Transverse Ising Hamiltonian 



N 

H = -rj^ erf - 2J JijViVj 

= H +V (17) 
The canonical partition function of TL reads 

Z = Tre~^ +v \ 

Now we apply the Trotter formula 

exp (Ai + A 2 ) = lim [expAi/M cxp A 2 /M] M , 
even when [A±, A 2 ] ^ 0. On application of this, Z reads 

Z = Y lim ( Si \ [exp (-PHo/M) exp {-(3V/M)] M \ Si ). (18) 



Here Sj represent the i-th spin configuration of the whole system, and the above summation runs over all such possible 
configurations denoted by i. Now we introduce M number of identity operators 



l= z 2\s i ,k){s i ,k\, fc = l,2,...M. 



in between the product of M exponentials in Z, and have 

M 



Z = ^lim^Tr JJ((Ti ifc ...o-jv,fc|exp ^ J ex P ) l^k+i-VN^+i), 
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and periodic boundary condition would imply <7jv+i, P = <?i, P - Now, 



M 



J|(o-i >fe ...crjv,fc|exp ( jjJ2 a * a 3 ) l°'i,fe+i- 'JV,fc+i) 



fc=i 



= cxp 



N M pj t - 



ij—l k—1 

where = ±1 are the eigenvalues of o z operator. Also, 



M 



JJ(CTi,fc...<Tjv,fc| exp 



k=i 



M 



|0'l 1 fe+l...(Tjv,fe+l} 



(19) 



= ( - sinh 



2/?r 



M 



exp 



JV M 



I Incoth (§)EE ^i.fc^.fc+i 



=i fe=i 



The last step follows because 



gacr 35 _ e -i(ia<T x ) 



cos (iacr x ) — i sin (iaa x ) = cosh (a) + cr 2 sinh (a), 



and therefore 



We°" \a'} = 



— sinh (2a) 



1/2 



exp [(<7<t'/2) Incoth (a)], 



since 



and 



(T 



T) = (l|e affa! |4>=cosh( a ) = 



(T|e aCT ||) = U \^\ T) = sinh (a) 



■ sinh (2a). coth (a) 



- sinh (2a) / coth (a) 



1/2 



1/2 



Thus the partition function reads 



Z = C™Tr a (-(m e ff[*}) ; C=isinh^ 



where the effective classical Hamiltonian is 

N M 



fe=l 



-jj<T lk <Jj k - ^ Incoth ^ ) on on. 



(20) 



(21) 



The Hamiltonian Tt e ff is a classical one, since the variables cafe's involved are merely the eigen- values of er z , and 
hence there is no non-commuting part in H e //. It may be noted from (21) that M should be at the order of h(3 (we 
have taken K = 1 in the calculation) for a meaningful comparison of the interaction in the Trotter direction with 
that in the original Hamiltonian (see Fig. 1.3). For T — > 0, M — > oo, and the Hamiltonian represents a system of 
spins in a (d+l)-dimensional lattice, which is one dimension higher than the original d-dimensional Hamiltonian, 
as is evident from the appearence of one extra label k for each spin variable (see Fig. 1.3). Thus corrosponding to 
each single quantum spin varible Oi in the original Hamiltonian we have an array of M number of classical replica 
spins o^ ■ This new (time-like) dimension along which these classical spins are spaced is known as Trotter dimension. 
From the explicit form of W e // ; w e see that in addition to the previous interaction (J) term (— JijVi&j), there 
is an additional nearest neighbour interaction (J') between the Trotter replicas corrosponding to the same original 
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FIG. 4: At the left is a portion of a one dimensional quantum Ising chain with random exchange interactions and at the right 
is a part of its Suzuki- Trotter equivalent classical lattice with randomness correlated in Trotter direction. 



spin, along the Trotter direction, given by the term . J^fcLi — lncoth (/3r/M)aikCTiK+i) ( a s shown in 

Fig. 1.3). For finite temperature, the optimal width of the lattice in the Trotter direction is finite and the critical 
behaviour remains c?-dimensional. 



The calculations, and consequently the effective Hamiltonian (21), is valid for any general interaction J^; ofcourse, 
r has been taken to be nonrandom. Fig. 1.4 describes a situation where Jij were nonrandom (we had Jij = J). 
For random Jy, where Jij were nonrandom (we had Jij = J). For random Jy, remain identical (</') wheras the 
spatial randomness in interactions for various Trotter slices get correlated as indicated in Fig. 1.4. Such equivalence 
of d-dimensional quantum system with a (d+ l)-dimensional classical model can also be seen from the renormalization 
group study of the quantum models (say, one-dimensional transverse Ising model and its equivalent critical behaviour 
of two-dimensional classical Ising system) as shown in Appendix B. 



VI. CLASSICAL SPIN GLASSES: A SUMMARY 



Spin glasses are magnetic systems with randomly competing (frustrated) interactions Q ■ Frustration is a situation 
where all of the spins present in the system cannot energetically satisfy every bond associated to them. Here the 
frustration arises due to competing (ferromagnetic and anti- ferromagnetic) quenched random interactions between 
the spins. As a result there arise huge barriers (O(N), N = system size) in the free-energy landscape of the system. In 
thermodynamic limit, height of such barriers occassionally go to infinity. These barrieres strongly separate different 
configurations of the system, so that once the system gets stuck in a deep valley in between two barriers, it practically 
gets trapped around that configuration for a macroscopically large time. Because of frustration, the ground state 
is largely degenerate; degeneracy being of the order of exp(iV). As discussed above, these different ground state 
configurations are often separated by O(N) barriers, so that once the system settles down in one of them, it cannot 
visit the others equally often in course of time, as predicted by the Boltzmann probability factor. The system thus 
becomes "nonergodic" and may be described by a nontrivial order parameter distribution 8] in the thermodynamic 
limit (unlike the unfrustrated cooperative systems, where the the distribution becomes trivially delta function-like). 
The spins in such a system thus get frozen in random orientations below a certain transition temperature. Although 
there is no long range magnetic order, i.e., the space average of spin moments vanishes, the time average of any spin 
is nonzero below the transition (spin-glass) temperature. This time average is treated as a measure of spin freezing 
or spin glass order parameter. 

Several spin glass models have been studied extensively using both analytic and computer simulation techniques. 
The Hamiltonian for such models can be written as 

H = -J2j ij af* z j (22) 

i<j 

where Sf = ±1,2, ...,N, denote the Ising spins, interacting with random quenched interactions Jij, which differs in 
various models. We will specifically consider three extensively studied models. 



(a) In Sherrington-Kirkpatrick (S-K) model Jij are long-ranged and are distributed with a Gaussian probability 
(centered around zero), as given by 

I N V' 2 f-NJ?,\ 

p < j «K^) exp (^j (23 ' 
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(b) In Edward- Anderson (EA) model, the Jy's are short-ranged (say, between the nearest neighbours only), but 
similarly distributed with Gaussian probability (23) 

(c) In another kind of model, the J^ 's are again short-ranged, but having a binary (±J) distribution with probability 
p: 

P(J t3 ) = P 8{J l3 - J) + (1 - p)S(J ij - J). (24) 

The disorder in the spin system being quenched, one has to perform configurational averaging (denoted by overhead 
bar) over hiZ, where Z{= Trexp— f3TL) is the partitation function of the system. To evaluate (InZ), one usually 
employs replica trick based on the representation InZ = limn^oK-^" ~ I)/"] - Now for classical Hamiltonian (with all 
commuting spin components), Z n = L|™ =1 Z a = ^Q2a=i H a ), where H a is the a-th replica of the Hamiltonian H. 
in equation (22) and Z a is the corrosponding partition function. The spin freezing can then be measured in terms of 
replica overlaps, and Edward- Anderson order parameter takes the form 

i N i N 

9 = ^E (summit (W>> 

i=l i=l 

where a and /3 corrosponds to different replicas. 



Extensive Monte Carlo studies, together with the analytical solutions for the mean field of S-K and EA models, 
have revealed the nature of spin glass transition. It appears that the lower critical dimension df for EA model, below 
which transition ceases to occur (with transition temperature T c becoming zero), is between 2 and 3: 2 < d\ < 3. 
Thu upper critical dimension d^, at and above which mean field results (e.g., those of S-K model) apply, appears 
to be 6: d c u = 6. Within these dimensions (df < d < d c u ), the spin glass transitions occur (for Hamiltonians with 
short-ranged interactions) and the transition behaviour can be characterized by various exponents. Although the 
linear susceptibility shows a cusp at the transition point, the nonlinear susceptibility XSG = 0-/N) ^ r g(r), where 
g(r) — (1/iV) Yli ({St 'Si+ r )) 2 ! ) diverges at the spin glass transition point : 

X sg~(T-T c )-^, g(r)~r-( d - 2 +^f(~y, £ ~ \T - T c \-"* (25) 

Here £ denotes the correlation length which determines the length scaling in the spin correlation function g(r) (/ in 
g(r) denotes the scaling function). Numerical simulation gives v c — 1.3 ± 0.1, 0.80± 0.15, 1/2 and 7 C = 2.9 ± 0.5, 
1.8± 0.4, 1 for d — 2,3 and 6 respectively for the values of exponents. One can define the characteristic relaxation 
time r through the time dependence of spin auto-correlation 

q(t) = (S*(t)S?(0))~t-*q(~y, r~$'~|T-T c |-"«/*«> (26) 

where x = (d— 2 + r; c )/2z c , and z c denotes the classical dynamical exponent. Numerical simulations give z c — 6.1±0.3 
and 4.8±0.4 in d = 3 and 4 dimensions respectively. Of course, such large values of z c (particularly in lower dimensions) 
also indicates the possibility of the failure of power law variation (26) of t with T — T c and rather suggests a Vogel- 
Fulcher like variation: r ~ exp [A/ (T — T c )]. In the ± J spin glasses (type (c) above), some exact results are known 
along the 'Nishimori Line' 0, and the nature of the phase transition there is precisely known. 



VII. QUANTUM SPIN GLASSES 



Quantum spin glasses |9l|-|l3j| have the interesting feature that the transition in randomly frustrated (competing) 
cooperatively interactimg systems can be driven both by thermal fluctuations or by quantum fluctuations. Quantum 
spin glasses can be of two types: vector spin glasses introduced by Bray and Moore (see Q), where of course quantum 
fluctuation cannot be tuned, or a classical spin glass perturbed by some tunable quantum fluctuations e.g., as induced 
by a non commutative transverse field 0,0- The amount of quantum fluctuation being tunable, this Transverse Ising 
Spin Glass (TISG) model is perhaps the simplest model in which the quantum effects in a random system can be and 
has been studied extensively and systematically^, ^| . Precise realization of TISG in LiHo x Yi_ x F4, with magnetic 
Holonium ion concentration around x = 0.167 [l2l |. has led to several important developments. 

The interesting in such quantum spin glass models is about the possibility of tunnelling through the (infinitely 
high) barriers of the free energy landscape in the classical spin glass models (e.g., S-K model) due to the quantum 
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fluctuations induced by the transverse field. In classcal system, the overriding of an infinitely high barrier is infinitely 
hard for thermal fluctuations at any finite temperature. But quantum fluctuation can make a system tunnel through 
such a barrier, if its width is infinitessimally small. The barrier widths are actually seen to decrease with system size 
indicating to an ergodic (replica symmetric) picture for the free-energy landscape. 



A. Models 



1. Sherrington- Kirkpatrick Model in a Transverse Field 
The sherrington-Kirkpatrick (S-K) model in presence of a non-commutating tunnelling field, given by the Hamiltonian 

ij i 

where the follows the Gaussian distribution 

/ N \ 1/2 f-NJ 2 \ 

P < J «> = G^) " P (^J (28) 

was first studied by Ishi and Yamamoto . 
Phase Diagram 

Several analytical studies have been made to obtain the phase diagram of the transverse Ising S-K model (giving 
in particular the zero-temperature critical field). The problem of S-K glass in transverse field becomes a nontrivial 
one due to the presence of noncommuting spin operators in the Hamiltonian. This leads to a dynamical frequency 
dependent (spin) self-interaction. 

(i) Mean field estimates : 

One can study an effective spin Hamiltonian for the above quantum many body system within the mean field frame- 
work. A systematic mean field theory for the above model was first carried out by Kopec (see e.g., Q), using the 
thermoficld dynamical approach and the short time approximation for the dynamical spin self-interaction. Before 
going into the discussion of this approach, we shall briefly review the replica-symmetric solution of the classical S-K 
model (r = 0) in a longitudinal field given by the Hamiltonian 

H = -E^'^M- /i E (r i ( 29 ) 

(ij) 

where follows the Gaussian distribution given by (56). Using the replica trick, one obtains for configuration 
averaged n- replicated partition function Z n , given by 

/CO 
P{Jii)dJu exp LflE J « E °te°ia +/3hJ2 < 

Performing the Gaussian integral, using Hubbard-Stratonovich transformation and finally using the method of steepest 
descent to evaluate integrals for thcrmodynamically large system, one obtains free energy per site /, given by 



-/3f = lim 



/3A 2 



U-iE^+^ nT Kex Pj L)j 



where L — (f3J) 2 J2 a p q a f3&a a i3 + ^ELi an( ^ lap is self-consistently given by the saddle point condition 
(df /dq a p) — 0. Cosidering the replica symmetric case (q a p = q), one finds 

(BA) 2 1 f°° ^ 

-pf =¥—L(l- q )+ dr e~-ln[2cosh{/3/i(r)}] 
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where r is the excess static noise arising from the random interaction Jy and the spin glass order parameter q is 
self-consistently given by 



1 f°° r 2 

q=—= dr e~~ tanh 2 {(3h(r)} 
V2n J-oo 



and h(r) = A^/qr + h can be interpreted as a local molecular field acting on a site. Different sites have different fields 
because of disorder, and the effective distribution of h{r) is Gaussian with mean and varience A 2 q. 

At this point we can introduce quantum effect through transverse field term —T J^. af (with longitudinal field h — 0). 
The effective single particle Hamiltonian in the transverse Ising quantum glass can be written as 

n s = -h z (r)a z -Ta x , 

where h z (r) 1 as mentioned earlier, is the effective field acting along the z direction arising due to nonzero value of the 
the spin glass order parameter. Treating h z (r) and a as classical vectors in pseudo-spin space, one can write the net 
effective field acting on each spin as 



ho(r) = h z (r)z ~ Tx; \h (r)\ = ^h z (rf+T 2 . 
One can now arrive at the mean field equation for the local magnetisation, given by 

m(r) = p(r) tanh [f3h Q (r)] ; p(r) = £-£tt - 

M r )l 

and consequently, the spin glass order parameter can be written as 



1 f 00 

t= / dr e- r2 / 2 tanh 2 {(3h (r)}p 2 (r) . 



The phase boundary can be found from the above expression by putting q — > (h z (r) — J ^fqr an d = r), when 
it gives 

= tanh [ -—] . (30) 

A \k B TJ v ' 

From above we get T c = J. Ishi and Yamamoto used the "reaction field" technique to construct "TAP" like equation 
for free energy of the Hamiltonian (27) and perturbatively expanded the free energy in powers of T upto the order T 2 
to obtain 



k B T c = A[l- 0.23(17 A) 2 ]. 



(ii) Monte Carlo Studies 



Several Monte Carlo studies have been performed f° r S-K spin glass in transverse field. Applying Suzuki- Trotter 

formulation (as discusseed earlier) of effective partition function, one can obtain the effective classical Hamiltonian in 
Mth Trotter approximation as 



Heff = ~M l^Zs J H a i^3k ~ 777, InCOth T7 



M ^ ^ J J 2(3 \M, 

i,j=l k=l ' v ' i=l k=l 



ik+1 



NM , 

In 



1 2/3T 
— smn 

2 M 



(31) 



where cr^ denotes the Ising spin defined on the lattice (£, fc), i being the position in the in the original S-K model 
and k denoting the position in the additional Trotter dimension. 



Ray et al [Hj took T << J and their results indeed indicate a sharp lowering of Tc(T). Such sharp fall of T^r) 
with large T is obtained in almost all theoretical studies of the phase diagram of the model. 
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2. Edward-Anderson Model in a Transverse Field 



The Hamiltonian for the Edward- Anderson spin glass in presence of transverse field is that given by (27), where 
the random interaction this time is restricted among the nearest neighbours and satisfies a Gaussian distribution with 
zero mean and variance J, as given by 




With r = 0, the above model represents the E-A model with order parameter q = {of) 2 — 1 (at T = 0). When the 
transverse field is introduced, q decreases, and at a critical value of the transverse field the order parameter vanishes. 
To study this quantum phase transition using quantum Monte Carlo techniques, one must remember that the ground 
state of a ci-dimensional quantum model is equivalent to the free enery of a classical model with one added dimension 
which is the imaginary time (Trotter) dimension. The effective classical Hamiltonian can be written as 

W = ^2 KijCik^jk — ^ ^ K(T ik Oik+l, (32) 
k ij k i 



with 



Kn = ^-x tf = rlncothf^ 



M 2 V M 



where are classical Ising spins and denotes the original d-dimensional lattice sites and k = 1, 2, M denotes 
a time slice. Although the equivalence between classical and the quantum model holds exactly in the limit M — ► oo, 
one can always make an optimum choice for M. The equivalent classical Hamiltonian has been studied using standard 
Monte Carlo technique. The numerical estimates of the phase diagram etc. are reviewed in details in [l3| 

B. Replica Symmetry in Quantum Spin Glasses 

The question of existence of replica-symmetric ground states in quantum spin glasses has been studied extensively in 
recent years. Replica symmetry restoration is a quantum phenomena arising due to the quantum tunnelling between 
the classically 'traped' states seperated by infinitely high (but infinitessimally narrow) barriers in the free energy 
surface, which is possible as the tunnelling probability is proportional to the barrier area, which remains finite. To 
investigate this aspect of quantum glasses, one has to study the overlap distribution function P(q) given by 



P(<z) = ^P*P^(g -?("')), (33) 
1,1' 

where P; is the Boltzman weight associated with the state I and q 11 is the overlap between the sates I and I' 

1 N 

4=1 

One can also define the overlap distribution in the following form (for a finite system of size N) 



P N (q) = (S(q-qW)), (35) 

where q( 12 ^ is the overlap between two sets of spins uj 1 ' and cj 2 \ with identical bond distribution but evolved with 
different dynamics, 



N 



4EW (36) 



Pn(q) — * P(q) m the thermodynamic limit. In quantum glass problem one can study similarly this overlap distridu- 
tion Pjv(g); and if the replica symmetric ground states exists, the above function must tend to a delta function in 
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thermodynamic limit. In para-phase, the the distribution will approach a delta function at q = for the infinite system. 

Ray, Chakrabarti and Chakrabarti jlOj . performed Monte Carlo simulations, mapping the d-dimensional transverse 
S-K spin glass Hamiltonian to an equivalent (d + l)-dimensional classical Hamiltonian and addressed the question of 
stability of the replica symmetric solution, with the choice of order parameter distribution function given by 



/ n m \ 

*M9) = <* U-j^EE'M >, (37) 

V i=i k=i / 

where, as mentioned earlier, subscripts (1) and (2) refer to the two identical samples but evolved through different 
Monte Carlo dynamics. It may be noted that a similar definition for q (involving overlaps in identical Trotter indices) 
was used by Guo et al |ll|. Lai and Goldschmidt performed Monte Carlo studies with larger system size (N < 100) 
and studied the order parameter distribution function 



PN(q) = (s{q- W Y ( all ) cT§j), (38) 

where the overlap is taken between different (arbitrarily chosen) Trotter indices k and k'; k ^ k'. Their studies 
indicate that -P/v(g) does not depend upon the choice of k and k' (Trotter symmetry). Rieger and Young (see Q) also 
defined q( 12 ' in similar way (q( 12 ^ = (l/NM) y^ r Ylw) a ik ° dfay There are striking differences between the results 
Lai and Goldschmidt obtained with the results of Ray et al yuj. For T « L c , P(q) is found to have (in 0]) an 
oscillatory dependence on q with a frequency linear in N (which is probably due to the formation of standing waves 
for identical Trotter overlaps). However, with increase in N, the amplitude of oscillation decreases and the magnitude 
of P(q = 0) decreases, indicating that P(q) might go over to a delta function in thermodynamic limit. The envelope 
of this distribution function appears to have an increasing P(q = 0) value as the system size is increased. Ray et 
al |10| argued that the whole spin glass phase is replica symmetric due to quantum tunnelling between the classical 
trap states. Lai and Goldschmidt on the other hand, do not find any oscillatory behaviour in P(q). In contrary they 
get a replica symmetry breaking (RSB) in the whole spin glass phase from the nature of P(q), which in this case, has 
a tail down to q = even as N increases. According to them their results are different from Ray et al [ljj because of 
different choices of the overlap function. Goldschmidt and Lai have also obtained replica symmetry breaking solution 
at first step RSB, and hence the phase diagram. 

Biittner and Usadel (see e.g., Chakrabarti et al Q), have shown that the replica symmetric solution is unstable 
for the effective classical Hamiltonian (58) and also estimated the order parameter and other thermodynamic 
quantities like susceptibility, internal energy and entropy by applying Parisi's replica symmetry breaking scheme to 
the above effetive Hamiltonian. Using static approximation, Thirumalai et al (see Q), found stable replica symmetric 
solution in a small region close to the spin glass freezing temperature near the phase boundary. But as mentioned 
earlier, in the region close to the critical line, quantum fluctuations are subdued by the thermal fluctuations. Thus 
the restoration of replica symmetry breaking, which is essentially a quantum effect, perhaps connot be prominent there. 

All these numerical studies are for equivalent classical Hamiltonian, obtained by applying the Suzuki- Trotter 
formalism to the original quantum Hamiltonian, where the interactions are anisotropic in the spatial and Trotter 
direction and the interaction in the Trotter direction becomes singular in the limit T — > 0. Obviously one cannot 
extrapolate the finite temperature results in zero temperature limit. The results of exact diagonalization of finite 
systems (N < 10) at T = itself do not indicate any qualitative difference in the behaviour of the (configuration 
average) mass gap A and the internal energy E g from that of a ferromagnetic transverse Ising case, indicating the 
possibility that the system might become "ergodic". On the other hand, the zero temperature distribution for the 
order parameter does not appear to go to delta function with increasing N as is clearly found for the corresponding 
ferromagnet (random long range interaction without competition). In this case the order parameter distribution P(q) 
is simply the number of ground state configuratons having the order parameter value as q. This perhaps indicate 
broken ergodicity for small values of T. The order parameter distribution also shows oscillations similar to that 
obtained by Ray et al [To| . 

Kim and Kim pTj have very recently investigated the S-K model in transverse field using imaginary time replica 
formalism, under static approximation. They have shown that the replica-symmetric quantum spin glass phase is 
stable in most of the area of the spin glass phase, as have been argued by Ray et al, in contrary to the results of Lai 
et al and Thirumalai et al (see e.g., Chakrabarti et al 3). 
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VIII. QUANTUM ANNEALING 
A. Multivariable Optimization and Simulated Annealing 

Multivarable optimization problems consists of finding the maximum or minimum values of a function (known as 
cost function) of very many independent variables. A given set of values for the whole set of independent variables 
defines a configuration. The value of the cost function depends on the configurations, and one has to find the 
optimum configuration that minimizes or maximizes the cost function. The explicit evaluation of the cost function 
for all possible configurations in this context, generally turns out to be absolutely impracticable for most systems. 

One can therefore start from an arbitrary state and go on changing the configuration following some stochastic 
rule, unless an extremum is reached. For example, in a minimization problem, one may start from an arbitrary 
configuration, change the configuration according to some stochastic rule, evaluate the cost function of the changed 
configuration, and then compare its value with that of the original configuration. If the new cost function is lower, the 
change is retained, i.e., the new configuration is accepted. Otherwise the change is not accepted. Such steps may be 
repeated for times unless a minimum is reached. But in most cases of multivariable optimisation problem, there are 
many local extrema in the cost function landscape, and one cannot be sure that the extremum that has been reached 
is the global one. Kirkpatrick et al 15] proposed a very ingenious physical solution to this mathematical problem, 
now known by the name simulated annealing. The basic underlying principle of simulated annealing as follows. It is 
known that an ergodic physical system, at any finite temperature resides in the global minimum of its free energy. 
The minimum of the free energy is a thermodynamic macro-state corresponding to a maximum number of accessible 
microcsopic configuration. Hence at thermal equilibrium an ergodic system explores its configuration space randomly 
with equal apriori probability of visiting any configuration, and consequently is found most of the time at one or other 
of the configurations that corrosponds to the free energy minimum (since the number of configurations corresponding 
to such minimum is overwhelmingly large compared to that of any other macro-state). Now if the system starts from 
an arbitrary macrostate (not the minimum of free energy) then due to thermal fluctuation it reaches the free energy 
minimum within some time r known as the thermal relaxation time of the system. 

For an ergodic system (away from critical point) this relaxation time increases linearly with system size (which 
is logerithmically smaller a number compared to the corresponding number of all possible configurations). Hence if 
one follows the random dynamics of the thermal relaxation of a system, then he will be able to reach the minimum 
of cost function (zero temperatur free energy) in a substantially smaller time. What one needs to do is to view the 
cost function E as the internal energy of some system and start from an arbitrary configuration. Then one changes 
the configuration according to somestochastic rule, just as before. Now if the energy is lowered by the change, the 
change is accepted, but if it is not, the change is not thrown away with certainity. Instead it is accepted with a 
probability equal to the Boltzmann factor e - AE / kBT : where AE = Euft eT change) — -E(bcforo change) (since this is the 
way how systems relax thermally to their free energy minimum). Temperature T here is an artifically introduced 
parameter which has a high value initially, and is reduced slowly as time goes on, finally tending towards zero. At 
zero temperature the free energy is nothing but the internal energy of the system, and thus at the end of the final 
stage of annealing the system can be expected to be found, with a very high probability, in a configuration that 
minimizes the internal energy (cost function). 

However this simulated annealing technique can suffer severe set back when the system is "nonergodic" , like the 
spin glasses we discussed earlier. In such cases configurations corrosponding to minimum of the cost function are 
separated by 0{N) sized barriers, and at any finite temperature thermal fluctuations will take practically infinite 
time to relax the system to the global minimum crossing these barriers in thermodynamic limit N — ► oo. 

B. Ergodicity of Quantum Spin Glasses and Quantum Annealing 

The non-ergodicity problem makes the search of the ground state of a classical spin glass a computationally hard 
problem (no algorithm bounded by some polynomial in system size exists for such NP-hard problems). The problems 
of simulated annealing of spin glass- like systems can be overridden (atleast partially) by employing the method of 
quantum annealing |l6l Il7|. The basic idea is as follows: First the problem has to be mapped to a corrosponding 
physical problem, where the cost function is represented by some classical Hamiltonian (say TLq) of the form (22). 
Then a suitably chosen noncommuting quantum tunnelling term (say H'(t)) is to be added so that the Hamiltonian 
takes the form of (27). One can then solve the time dependent Schrodinger equation 



(39) 



15 




configurations > 



FIG. 5: Schematic indication of the advantage of quantum annealing over classical annealing. 



for the wave-function ip(t) of the entire system Ha + Ti'(t). The solution of the time dependent schrodinger equation 
approximately describes a tunnelling dynamics of the system between different eigenstates of 7io- Like thermal 
fluctuations in (classical) simulated annealing, the quantum (tunnelling) fluctuations owing to TL' in (39) help the 
system to come out of the local 'trap' states. If TC'(t) — > for t — > oo, the system eventually settles in one of the 
eigenstates of 7io! hopefully the ground state. The introduction of such a quantum tunnelling is supposed to make 
the infinitely high (but infinitessimally thin) barriers transparent to the system (see, e.g., Appendix C), and it can 
make transitions to different configurations trapped between such barriers, in course of annealing. In other words, 
it is expected that application of a quantum tunnelling term will make the free energy landscape ergodic, and the 
system will consequently be able to visit any configuration with finite probability. Finally the quantum tunnelling 
term is tuned to zero {Ti.'(t) — > 0) to get back the classical Hamiltonian. It may be noted that the success of quantum 
annealing is directely connected to the replica symmetry restoration in quantum spin glass [Tol ] , |14| due to tunnelling 
through barriers (see Fig. 1.5 and the discussion in the preceeding section). 

Here, the d-dimensional quantum Hamiltonian (27) (to be annealed) is mapped to the (d+ l)-dimensional effective 
Hamiltonian 



M / N n \ 

fc+l \ i,3 i=1 j 

where 

mt ( r \ 

J 1 = lntanh > 

2 \MT) 

is the nearest neighbour ferromagnetic coupling in Trotter direction, between the Trotter replicas of the same spin. In 
course of annealing, the temperature is kept constant at a low but nonzero value, and the tunnelling field T is tuned 
slowly from a high initial value to zero. The decrease in T amounts to the increase in J' (as casn be seen from above 
expression of J'). Initially at high T, J' is low, and each Trotter replica behaves almost like an independent classical 
spin system. The tunnelling field is then lowered in small steps. In each such step, the system is annealed in presence 
of the small temperature. Finally as T — > 0, J' — > oo, forcing all replicas to coinside at the end. As mentioned already, 
quantum annealing possibility directly rests on the replica symmetry restoration feature of quantum spin glasses, 
discussed in earlier section. In fact ergodicity in quantum spin glasses, as suggested in Fig. 1.5 was attributed by Ray 
et al to the "quantum fluctuations due to transverse field. Quantum tunnelling between the classical 'trap' states, 
separated by infinite (but narrow) barriers in the free-energy surface, is possible as quantum tunnelling probability is 
proportional to the barrier area which is finite." (Sec. V, Ref. [lfj|). 
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C. Quantum Annealing in Kinetically Constrained Systems 

It is largely believed that apart from the complexity associated to the non-trivial ground state structure of a glassy 
system, the occurance of certain kinetic constraints (blockings) during relaxation also contributes substantially to the 
slowness of its low temperature dynamics. The Kinetic constraints or blockings can be viewed as infinitely high energy 
barriers appearing in the relaxation path of the system. In order to relax to the minimum of the free energy, the 
system has to jump over these high barriers thermally, which they fail to do at any finte temperature. However if such 
barriers are infinitessimally narrow, then the system might be able to tunnel through them quantum mechanically 
if sufficient quantum fluctuation T is present in the system. Thus if one tries to anneal such a system down to its 
ground state starting from an arbitrary state, then quantum annealing might turn out to be much superior to the 
thermal annealing (see e.g., |l8Lll9|). 

We have studied [23] the annealing of a kinetically constrained Ising spin chain of N spins, starting from a dis- 
ordered state (with negligible initial magnetization), to its (external field induced) fully ordered ground state. At 
any finite temperature T (in the classical model) the system takes an exponentially long time to relax to the ordered 
state because of the kinetic constraints, which act like an infinite potential barrier, depending on the neighbouring 
spin configurations. Quantum mechanically, this infinite barrier is taken to be penetrable, ie with finite tunnelling 
probability, depending on the barrier height x an d width a (a — > faster than x~ 2 )- The introduced noise, required 
for the annealing, is reduced following an exponential schedule in both the cases: T — Toe^ l ^ c , T — Toe^ t / TQ , with 
Tq To. For our simulation for the quantum case, we have taken the tunnelling probabilities P (for cases I- IV) and 
employed them in a semi-classical fashion for the one dimensional spin chain considered. We observe that for similar 
achievement in final order (m/ ~ 0.92 starting from m. l = 1(T 3 ), r c ~ 10 3 tq for N = 5 x 10 4 . For even larger order 
(nif ~ 1), quantum annealing works even better (tq ~ 10 3 tq, for the same value of N). These comparison are for 
g = 10 2 and x — 10 3 for the constraint barriers [20]. 

In this picture, we considered the collective dynamics of a many particle system, where each one is confined in 
a (field) induced asymmetric double well potential for which we considered only the low lying two states (the wave 
packet localized in one well or the other), representing the two states (up and down) of an Ising spin discussed above. 
The tunnelling of the wave packet from one well to the other was taken into account by employing a scattering picture 
and we used the tunnelling probabilities as the flip probabilities for the quantum Ising spins. As such, the reported 
simulation for the one dimensional quantum East model is a semiclassical one. It may be noted however that, because 
of the absence of inter-spin interaction, the dimensionality actually plays no role in this model except for the fact that 
the kinetic constraints on any spin depend only on the left nearest neighbour (directedness in one dimension). Hence 
the semiclassical one dimensional simulation, instead of a proper quantum Monte Carlo simulation (equivalent to a 
higher dimensional classical one |^), is quite appropriate here. 



IX. SUMMARY AND DISCUSSIONS 



We have introduced the transverse Ising model for discussing the order-disorder transition (at zero temperature) 
driven by quantum fluctuations. Mean field theories are discussed next in sections 1.3 and 1.4. Application to BCS 
superconductivity theory is discussed in appendix A. Renormalization group technique for study of critical behaviour 
in such quantum systems is discussed in appendix B (for a chain). Next we have discussed the Suzuki- Trotter mapping 
of the d-dimensional quantum system to d + 1 dimensional classical system (in section 1.6). We introduce then the 
transverse Ising spin glass models, namely, the S-K model in transverse field and the E-A model in transverse field 
(Sec. 1.8.2). The existing studies on their phase diagrams are discussed briefly. We then discuss about the problem 
of replica symmetry restoration in quantum spin glasses (in Sec. 1.8.3). The application of the quantum annealing 
technique to capture the near-global minima of NP-hard problems is then discussed, and the effectiveness of quantum 
tunnelling over the thermal barrier hopping is discussed (Sec. 1.9). 

It may be noted in this connection that some recent attempts have been made to apply similar annealing, induced 
by quantum fluctuations, to the optimization problems like the travelling salesman problem, image restoration, etc. 
hahj. Like the near-global minima in free energy landscape of such spin glasses, the barriers are often globally 
contributed and these barrier heights grow as the system size grows (unlike the locally optimized configurations and 
the barriers between them). Classically, the system becomes nonergodic due to these macroscopically high barriers 
(NP-hard to reach the ground state), as thermal fluctuations have to wait until they can scale such macroscopically 
high barriers. Quantum tunnelling does not necessarily look for barrier height |Toj to overcome them (by tunnelling 
through; see appendix C, see Fig. 1.5) and helps restoring replica symmetry as well as annealing [Trl Il7j. 
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X. APPENDIX 

Appendix A: 

MFT of TIM and BCS Theory of Superconductivity 

The phonon mediated effective attractive interaction between electrons give rise to a cooperative quantum Hamil- 
tonian. Although the quantum phase transition in such a system is not physical or meaningful, the finite tempera- 
ture superconducting phase transition can be studied easily following the mean field theory discussed here (using a 
pseudo-spin mapping The relevant part of the Hamiltonian of electrons that take part in superconductivity has 
the following form 

Here the suffix fc indicates a state with momentum fc and spin up, while (— k) indicates a state with momentum — fc 
and spin down and V is a positive constant that models the attractive coupling between electrons through phonons. 
We will solve this equation following spin-analog method Q . Here we are considering low- lying states containing pair 
of electrons (fc, —As). For a given fc, there are two possible states that come into consideration: either the pair exists, 
or it does not. Thus we enter into a spin-like two-state picture as follows. 

Let us introduce the number operator hk — c[.c?c . This reduces the Hamiltonian (Al) to 

Ti-red = ~y^efc(l - h k - n_ fe ) - V ^ 4' c -fc' c -fc c fc- (^ 2 ) 

fc kk' 

Here we have introduced a term — ^2 k e k with the choice ^2 k e& = in mind, for all fc's (basically these sums are 
over the states within energy ±lud about the fermi level, where is the Debye energy) that partictpates in pair 
formation. As stated earlier, if n k denotes the number of electrons in fc-state, then we are considering only a subspace 
of states defined by n k = where either the both of the states in the pair (fc,— fc) are occupied, or both are empyt 
. Now if we denote by Ilfel-fe) a (fc,— fc) pair-occupied state and by |0fe0_fc) an unoccupied one, then 

(1 - h k - n_fc)|lfel_fc) = (1 - 1 - l)|lfcl- fc ) = — ]l fe l_fc>, 

and 

(1 - h k - n_ fc )|0 fc 0_fc) = (1 - - 0)|0 fc 0_ fc ) = |0 fe 0_fc) 
Thus we switch over to our good old pseudo-spin picture through the following corrospondences 

|i*i-fc) I 
|o fc o_ fe ) «*■ I T) fc , 

and (1 — n k — ri-k) (^3) 

Since 

4 c -fcl T>fc = I l)k, c\ c -k\ i>fe = & c -kC k \ l)k = | T)fc, c_ fcCfc | T)fe - 0, 
we immediately identify its correspondence with raising and lowering operator a + / o~ : 

and 

a = { ) 
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and therefore 

c l c -k = \°k> c -kCk = (A4) 
Hence in terms of these spin operators, Hamiltonian (A2) takes the form 

k kk' 

Since the term X^fcfc' ( (T fc' cr fc — a k' a k) vanishes due to symmetric summing done over k and k', the Hamiltonian finally 
reduces to 

n = - E ^k \v E(^ + ( A6 ) 

k kw 

To obtain the energy spectrum of the pseudo-spin BCS Hamiltonian (A6) we apply now the mean field theory 
developed in earlier section. 

Weiss' Mean Field Solution 



Just as we did in case of TIM (see sec 1.3), here also we introduce an average effective field hfc for each pseudo-spin 
a k as 



2 

k> 

and conseqently the Hamiltonian (A6) takes the form 

k 

Here for each k there is an independent spin a k which interacts only with some effective field h^, and our system is a 
collection of such mutually non- interacting spins under a field . 

Now if redefine our x-axis along the projection of hk on the x-y plane for each k, then with all (of,) = we get 

tan^ = M = iZ^Ml , m 

K £ fe 

where 6 k is the angle between z-axis and hfe. 
Excitation spectra at T = 

Since at T = (a x ) = 1, 

= \a\ sm6 k , = sm9 k > (A8) 

Thus from (A7) we get 

tan0 fe = («/2e fc )y^sin 9 k > 

fe' 

Now let us define 



A= V^sin6» fe ,, 



2 

fe' 



so that tan#fc = A/e k . Then simple trigonometry gives 



an^-^j ; cosO k = -0=^ (A9) 
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Substituting for sin6>fc< into the above equation we get 

A 



2 ^-wt^; 

Assuming the spectrum to be nearly continuous, we replace the summation by an integral and note that V is attractive 
for energy within ±ud on both sides of fermi level; lod being of the order of Debye energy, then the last equation 
becomes 



l = -V PF j_^ 7 == = V PF ^-\. D /A). 



Here p F is the density of states at fermi level. Thus 

A= . J D =2cu D e- 1 / v ^ 1 ( if PF V«l) (A10) 
smh(l/Vp F ) 

We see that A is positive if V is positive. To interprate the parameter A, one may notice that at first approximation, 
the excitation spectrum is obtained as the energy £ k to reverse a pseudo- spin in the field h fc , i.e., 

£ k = 2\h k \ = 2(4 + A 2 ) 1/2 . (All) 

From this expression we clearly see that the minimum excitation energy is 2A, i.e. A gives the energy gap in the 
excitation spectrum. 

Estimating transition temperature T c 



To find the critical temperature for BCS transition, we just extend here the non-zero temperature version of mean 
field theory done for Ising case. We should have (unlike that in (All), where (a k ) = 1) for T = 0: 

(<4) = tanh (/?| h fe |) . (A12) 

Equation (23) accordingly modifies to 

tanfl fe = (^Wtanh(/?|hfc,|)sin0fc, = (A13) 

^ k ' k' £fe 

where A(T) = y J^k' tanh (^nf^j sin^/. From equation (All) we have 

|h fe |= [e£ + A 2 (T)]. 

The BCS transition is characterized by the vanishing of the gap A, since without such a gap in the spectrum, infinite 
conductance would not be possible except at T = 0. Hence, as T — > T c , A — > 0, i.e., using (All), |hfc| = e k and 
putting this and relations like (A9) in (A13), we get 

Above relation is correct if we consider an excited pair as a single entity. However, if we extend our picture to incor- 
porate single particles excited symmetrically in momentum space, then we double the number of possible excitations, 
thereby doubling the overall entropy. This is exactly equivalent to a doubling of the temperature in free energy. The 
energy contribution to the free energy, however, remains unaltered, since two single particle excitations of same |k| 
have same energy as that of a pair of equal |k|. Hence we replace T c by 2T Cl and in the continuum limit, get 

UD de , / e \ „ r°l 2T - tanhx , 
tanh — — = 2 / ax, 



V PF J- UD e \2T C 
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with (x = e/2T c ) This is the equation from which we obtain T c on integration. If T c << ujd, then we may approximate 
tanhx rs 1, for x > 1, and for a; << 1, we set tanhx ~ x. This readily reduces the integral to the value l+log(u>£>/2T c ), 
from which we have 

T c = (e/2)cj D e-V v ?». 

Grphical integration gives a closer result 

T c = l.Uu D e- 1/vpF . (A1S) 

Comparing equations (A10) and (A15) we get the approximate relationship 

2A ~ 3.5T C . (A16) 

This result is quite consistent with the exprimental values for a number of materials. For example, the value of 2A/T C 
are 3.5, 3.4, 4.1, 3.3 for Sn, Al, Pb, and Cd superconductors respectively. 

Appendix B: 

Real Space Renormalization for Transverse Ising Chain 

Here the basic idea of real space block renormalization Q, Q is illustrated by applying it on an Ising chain in 
transverse field. Taking the cooperative interaction along x-axis, and the transverse field along z-axis, the Hamiltonian 
reads 



N N-l 

~i u i+l 



h = -r£>?- J ^2 of of 



Here 



and 



= H B +H IB (say). (SI) 

N/b b 6-1 

H B = ^H P ; Hp = - ^ T(T Ip ~^2 J a i, P a i+i,p ( m ) 

p—l i—l i—1 

N/(b-l) 

Hib = ; T~Lp,p+i = _ ^ cr 6,p cr i,p+i- (-S3) 

P =i 

The above rearrangement of the Hamiltonian recasts the picture of N spins with nearest-neighbour interaction into 
one in which there are N/(b— 1) blocks, each consisting of b number of spins. The part TIb represents the interaction 
between the spins within the blocks, while TIib represents the interactions between the blocks through their terminal 
spins(see Fig. 1.6). 

Here we will consider b = 2, as shown in the figure. Now Ti. p has got 4 eigen-states, and one can express them in 
terms of the linear superposition of the eigen-states of af ® erf p ; namely, 

ITT), III), III), and ||T). 
Considering the orthonormality of the eigen-states, one may easily see that the eigenstates of Tt p can be expressed as 

10) = —L=(\rt) + a\ II)) 
Vl + a 

|i) = ^=(IU> + UT» 
|2) = ^=(|U>-UT» 
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°1,p °2,p °1,p+1 a 2,p+1 



• • 



p-th Block (p+1)th Block- 

H p,p+1 



FIG. 6: The linear chain is broken up into blocks of size 6 (= 2 here) and the Hamiltonian (Bl) can be written as the sum of 
block Hamiltonians TL V and inter-block Hamiltonians H PlP +i. The Hamiltonian Tip is diagonalized exactly and the lowest lying 
two states are identified as the renormalized spin states in terms of which the inter-block Hamiltonian is rewritten to get the 
RG recursion relation. 



Here a is a coefficient required to be chosen properly, so that |0) and |3) are eigenstates of H p . One gets, 

1 



H P \0) 



vT 



H)+a| II) 



[-IV? + <7f)-J(<7faf 

1 



1 



TT) + «I ID) 



r(2| H)-2a| ||))- J(\ii)) + a\ tt>] 



= -(2r + Ja) 



Thus |0) to be an eigenstate of H p , one must have 



1 



=> 



or, 



2T - J /a 
' 2T + J a 

J a 2 -4Ta-J = 

±V4r 2 + J 2 - 2T 



a = 



J 



To minimize the energy, we have to choose, 



One can now see, applying 7i p on its eigen-states, 



V4T 2 + J 2 - 2T 
J ' 



n p \o) 


= Eo\0), 


E 


H p \l) 


= E,\l), 


Ei 


n p \2) 


= E 2 \2), 


E 2 


n p \s) 


= E 3 \3), 


E3 



(B5) 



(B6) 
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Now we define our new renormalized spin variables er"s, each replacing a block in the original Hamiltonian. We retain 
only the two lowest lying states |0) and |1) of a block and define corresponding o'^ to have them as its two eigenstates, 
| T)' = |0) and | |> = |1). We also define 

, x af®l + l®a% 



where X is the 2x2 identity matrix. Now since 



we take our renormalized J to be 



J' = J 



1 + a 



V 2 (l + « 2 )' 

(i + q) 2 

2(1 + a 2 )' 



(B7) 



and since the energy gap between |0) and |1) must be equal to 2T' (This gap was 2r in the unrenormalized states), 
we set 



where a — V4A 2 + 1 — 2A., defining the relevant variable A = T/J. 
The fixed points of the recurrence relation (rewritten in terms of A) are 

A* = 
A* — > oo 
and A* ~ 1.277. 

Now if correlation length goes as 

c~(A-A c r, 

in the original system, then in the renormalized system we should have 

£' ~ (A' - x c f 



(B9) 



(B8) 



-> 



A' - A c 

X-Xr 



=> 



-X/u 



dX' 
~dX 



A=A C = A* 



(BIO) 



Now since the actual physical correlation length should remain same as we renormalize, £' (correlation length in the 
renormalized length scale) must be smaller by the factor b (that scales the length), than £ (correlation length in 
original scale), i.e., = b, or, 



h -x/u = 

dX 



O (say), 



A=A C =A* 



hence, 



TnT 



A=A* 



TnT 



1.47, (for b = 2), 



(BU) 



compared to the exact value v = 1 for (d+ 1 =) 2 dimensional classical Ising system. Similarly E g ~ ui ~ (time) 1 
£~ z ; z = 1. But for 6 = 2, we donot get z = 1. Instead, A'/A ~ gives z ~ 0.55. Energy gap 



A(A)-|A c -A| s -r z ^|A c -Ar z (B12) 
Hence s — vz — 0.55 x 1.47 ~ 0.81 (compared to the exact result s = 1). Results improve rapidly for large b values 
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1 



V = - 

FIG. 7: Quantum tunnelling thruogh asymmetric barrier. 



Appendix C: 

Tunnelling Through Asymmetric Barrier 



Let us consider an asymmetric potential energy barrier in one dimension, (as shown in Fig. F7). It is essentially a 
rectangular barrier of height x an d width a between two different energy levels with a potential difference h between 
them. The potential energy V is zero at the left of the barrier (x < 0), and it is h (h may be negative as well) at the 
right of the barrier (x > a). If a beam of free particles of mass m with kinetic energy T incidents on the barrier from 
the left, then one can calculate the probability for a particle in the beam to get transmitted through (or reflected by) 
the barrier by solving the time-independent Schrodinger equation (with a time- independent V). The transmission 
coefficient T (defined below) describes the probability of transmission for a single particle, as well as the average 
transmission of the incident beam. 

The incident wave function (x) , the intermediate wave function tp2 {%) and the transmitted wave function ip3 (x) 
then takes the form 



■01 (x) 


= Ae~ lklX 




= Be k2X + 




= De lk3X , 



x < 0, 

— k 2 x 

) 

x > a 



< x < a 



where, 



X 



and 



^3 



r-h, 



setting 2m/h — 1. Here A and D are the amplitudes of the incident and the transmitted wave respectively. At this 
point one may note that for T < h transmission is trivially zero. Hence we consider the case for T > h i.e., for real k^. 
In that case, applying the condition of continuiety of the wave function and its space derivatives at the boundaries, 
one obtains the relation (cf. pH*) 

A = \-De lk3a [{\ + k 3 /ki) cosh^a + i(n/ki — k 3 /n) sinhfta], 



r. We now consider the limit of very high but narrow barrier, such that \ ~ > °°? a — > 0, with 

X, and of course k is real. Since T > 0, fci is also real. Hence 



where k 2 = —k\ = \ 

g = x a finite. We also assume that \ ^> T, so that k 2 
under this condition the transmission coefficient defined as T = 

4fc 3 /fci 



\D\ 2 k s /\A\ 2 h is given by (cf. M) 



T 



1 + fa ) cosh^ (/so) 



fei 



- ^ ) sinh 2 (m) 
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In the limit of high but narrow barrier specified above, one has ko < 1. Hence neglecting terms quadratic or of higher 
order in na and linear in one gets 



T 



4fc 3 /fci 



(1 + 17) M 2 
^ 4yr(r - h) 
(VT+ v / r^) 2 + .g 2 ' 

putting k\ = Vr, ^3 = yT — h and K 2 a w = g. The transmission coefficient T is thus finite even when the barrier 
height x diverges keeping g = x a finite. 
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